REPORT No. 752 



ON A NEW METHOD FOR CALCULATING THE POTENTIAL FLOW PAST A 

BODY OF REVOLUTION 

By Cael Kaplan 



SUMMARY 

A new method is presented for obtaining the velocity potential 
of the flow about a solid of revolution moving uniformly in the 
direction of its axis of symmetry in a fluid otherwise at rest. 
This method is based essentially on the fact that the. form of the 
differential equation for the velocity potential is invariant with 
regard to conformal transformations in a meridian plane. By 
means of the conformal transformation of the meridian profile 
into a circle a system of orthogonal curvilinear coordinates is 
obtained, the main feature of which is that one of the coordinate 
lines is the meridian profile itself. The use of this type of 
coordinate system yields a simple expression of the boundary 
condition at the surface of the solid and leads to a rational 
process of iteration for the solution of the differential equation 
for the velocity potential. It is shown that the velocity potential 
for an arbitrary body of revolution may be expressed in terms 
of universal functions which, although not normal, are obtain- 
able by means of simple quadratures. 

The general results are applied to a body of revolution ob- 
tained by revolving a symmetrical Joukowski profile about its 
axis of symmetry. A numerical example -further serves to 
illustrate the theory. 

INTRODUCTION 

The simplest case of a three-dimensional fluid motion 
occurs when a hody of revolution moves with a constant 
velocity in the direction of its axis of symmetry. In this 
case, the. motion is the same in any plane passing through 
the axis of symmetry and, in this respect, presents some 
analogy with a two-dimensional motion. Thus, 'a stream 
function is defined hy means of the equation of continuity; 
and the condition for irrotational motion yields a velocity 
potential. The stream function and the velocity potential, 
however, are not interchangeable in three-dimensional flows 
in the same way as are the corresponding quantities in two- 
dimensional irrotational motions. The reason for this 
difference is that, although the differential equation for 
the velocity potential is Laplace's equation, the equation 
for the stream function is not Laplace's and therefore the 
two functions cannot be combined to give an analytic 
function of a single complex variable. It follows that the 
elegant and powerful methods of the complex variable 
are not obviously applicable and the calculation of the 
potential flow past a body of revolution has, of necessity, 
developed along other lines. 



The method of calculating the flow past a body of revolu- 
tion most often referred to was suggested by Rankine and 
developed by von Karman (reference 1) and others. For 
axial flow, the axis of the body is covered by a continuous 
distribution of sources and sinks in such a way that the closed 
stream surface found by the superposition of the flow in- 
duced by the sources and sinks on the parallel flow coincides 
with the surface of the solid. For transverse flow, the axis 
of the body is covered by a continuous distribution of dou- 
blets. The superposition of the flow due to the doublets on the 
parallel flow then yields the surface of the solid as a stream 
surface. Both problems lead to integral equations that von 
Kaim&n solved by an approximate method. As von Kar- 
man pointed out, however, the exact replacement of the body 
by a distribution of singularities along the axis of symmetry 
is possible only when the analytical continuation of the po- 
tential function, free from singularities in the space outside 
the body, can be extended to the axis of symmetry without 
encountering singular points. (See reference 1, p. 27.) In- 
asmuch as the question of when this analytical continuation 
is or is not possible has never been answered, numerical 
calculations may lead to incorrect results, particularly in 
the case of a body with a rather blunt nose. 

In reference 2 an attempt was made to calculate the po- 
tential flow past a body of revolution according to the 
methods of potential theory. In that paper, Laplace's equa- 
tion for the velocity potential is expressed in terms of elliptic- 
cylindrical coordinates and is solved in conjunction with 
the appropriate boundary conditions for axial and trans- 
verse flows. In the case of an ellipsoid or a hyperboloid of 
revolution the solution is obtained in a closed form, both 
shapes being members of the family of orthogonal coordi- 
nate surfaces belonging to a system of elliptic-cyhndrical 
coordinates. For any other solid of revolution, however, 
the method leads to two sets of linear equations, each set 
having an infinite number of equations and an infinite num- 
ber of unknown coefficients for the detennination of the 
velocity potentials for the axial and transverse flows. 

In the present paper, a new method is presented for cal- 
culating the potential flow past an arbitrary body of revolu- 
tion. Only the case of axial flow is discussed but the method 
is equally applicable to the case of transverse flow. The 
method is based on the discovery that, by the proper choice 
of the system of coordinates to be used for a given body 
of revolution, the solution of the potential-flow problem can 
be obtained by means of quadratures. Coordinate systems 
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of this nature already exist in the literature. For example, 
if the solid is an ellipsoid or a hyperboloid of revolution, the 
coordinate system used is an elliptic-cylindrical one. These 
coordinate systems possess in common the property that one 
of their coordinate surfaces is the boundary of the solid 
itself. It will be shown in this paper that, by means of the 
conformal transformation of the meridian profile into a 
circle, a system of orthogonal curvilinear coordinates can be 
denned such that one of the coordinate fines is the meridian 
profile of the solid of revolution. Furthermore, it will be 
shown that, by the use of this type of coordinate system, 
the potential-flow problem for an arbitrary body of revolu- 
tion can be solved by means of elementary integrations. 

In comparison of the method of this paper with the 
method of reference 2 it is to be recalled that in reference 2 
Laplace's equation for the velocity potential was expressed 
in elliptic-cylindrical coordinates. The general solution in- 
volves normal functions: namely, Legendre functions of the 
first and second kinds. The boundary condition, however, 
for- a$. arbitrary body of revolution involves the two inde- 
pendent variables of the problem and yields, as mentioned 
before, an infinite set of linear equations with an infinite 
number of unknown coefficients. On the other hand, the 
method of this paper utilizes a different coordinate system 
for each body of revolution. The main feature of these 
conformed orthogonal curvilinear coordinates is that one of 
the coordinate lines coincides with the meridian profile of 
the body of revolution. This fact leads to an expression 
of the boundary condition at the surface involving only one 
independent variable. Although the general solution of the 
differential equation for the velocity potential has not been 
obtained in terms of normal functions, a method of iteration 
has been devised that involves only elementary integrations. 
This solution, satisfying the boundary conditions, can then 
be expressed in terms of universal functions. Although 
these functions are not normal, they need he determined 
but once. 

MATHEMATICAL DEVELOPMENTS 

EQUATION FOR THE VELOCITY POTENTIAL 

The axis of symmetry of the body is denoted by x and the 
position of a point in a meridian plane is fixed by the Car- 
tesian coordinates (z, a) (fig. 1). Then if q x and q a are, 

Y 




Fioobe 1.— The meridian plane XOa 

respectively, the components of the fluid velocity* in the 
directions of the a; and a axes, the equation of continuity is 



obtained by equating to zero the flow out of the annular 
space obtained by revolving a small rectangle dx d& around 
the axis of symmetry. Thus, the total flow outward in the 

direction of the x axis is ^p,iraqJ[o)dx and in the direction 

of the gj axis is ^(2ir©gw£i;)cZ©. The equation of continuity 
is therefore 

l|(aff.)+|-(«2»)=0 (1) 

Since the flow is symmetrical about the x axis, the vor- 
ticity is 

dga Pgr 

Ox 0©' 

and if, further, the motion is irrotational, then 

bx do 

A velocity potential <£ can therefore be so defined that 

b<j>) 

and the equation of continuity (1) becomes 



(2) 



(3) 



Consider now the conjugate complex variables z=x-\-ia 
and z=x—ia. Then, symbolically, 

OZ 02 Otf 

*dz~dx ±% ba 
and it can be easily verified that 



dx 



(^)+s(»i)=^l(4§) 



Thus, the vanishing of the real part of ^ (°^f) 13 equivalent 

to equation (3) . 
Consider further the conformal transformation 
z=/(f) where £=£+171 

Then, 

Sinec z, s, or f, f may be looked upon as independent 
variables, it follows that the right-hand side of equation (4) is 



R.PA 



d^dfb 
dz dz 
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Now, the product ^ ^= is a real quantity. The vanishing of 

the real part of^^o^=) is therefore equivalent to equation (3). 
Hence, 

where ra is now a function of (£, »j) obtained from the con- 
formal transformation z=f(£). 1 

Equation (5) rather than equation (3) will be considered 
to be the fundamental differential equation governing axi- 
symmetrical motions of a perfect incompressible fluid. 
It is more general than equation (3) in the sense that the 
independent variables (£, rj) denote any- set of orthogonal 
curvilinear coordinates obtained by means of a conformal 
transformation z=j{£). In addition, its form is invariant 
with regard to conformal transformations and, in this sense, 
does not complicate the original equation (3) . In the follow- 
ing section it will be shown that, by means of the conformal 
transformation of the meridian profile of an arbitrary body 
of revolution into a circle, an orthogonal curvilinear system 
of coordinates (£, ij) can be so defined that the coordinate 
lino 77=0 is the meridian profile itself. In subsequent sec- 
tions it will be seen that the use of this type of coordinate 
system leads to a simple expression of the boundary condi- 
tion and, consequently, to an iteration process for the solu- 
tion of equation (5) for an arbitrary body of revolution. 

CONFORMAL TRANSFORMATION AND ORTHOGONAL COORDINATES 

It is well known that a unique conformal transformation 
exists which maps the region external to a given boundary 
in the z plane into the region external to a circle in the Z 
plane with its center at the origin, such that the regions at 
infinity of the two planes correspond. The function repre- 
senting this conformal transformation can be developed in 
the region external to the circle in a convergent series of the 
type 

2 =Z+ Cl ... (6) 

where the coefficients c u a lt a*, a*, . . . are, in general, com- 
plex quantities. In this paper, the meridian profile is sym- 
metrical with respect to the axis of revolution and these 
coefficients, therefore, are all real. The constants a u a 2 , 

1 It has been pointed out by Dr. Theodore Theodorsen that equation (5) can bo obtained 
directly, as follows: 
Tho vector form of equation (3) In a meridian plane Is simply 

dlv (n)irad<4)~n 

If a s et of orthogonal curvilinear coordinates { and q Is Introduced, (he expression for 
dlv (m grad <j>) becomes 

shore the eloments of length along the f-varlablo and ^-variable coordinate lines are, respec- 
tively, Ai df and ij dij. 

If, bowover, the transformation from the rectangular Cartesian coordinates to the 
orthogonal curvilinear coordinates (f,ij) is also conformal, then Ai~Ai; because for a conformal 
transformation tho magnification at any point In the plane is independent of the direction 
through the point. Thus, with Ai~A>, tho foregoing expression placed equal to zero yields 
equation (6). 



as. . . . depend only on the shape of the meridian profile in 
the s plane. 

The method of Theodorsen and Garrick, described in 
reference 3, is particularly well adapted for the purpose of 
determining transformation (6). It is shown in reference 3 
that, with the proper choice of axes and origin, the Joukowski 
transformation 

z=Z'+§ (7) 

maps a closed boundary in the z plane into a nearly circular 
curve in the Z' plane. The mapping of this nearly circular 
boundary, with its center at the origin of the Z plane, is. 
then completed by means of the transformation 

Z'=Zef™ (8) 

where 

n-l A 

On elimination of Z' between equations (7) and (8), equa- 
tion (6) follows, where 

<h=C2 +|c 1 2 +a* 

a?=c i +CiC z +^i s <h+^*+^S+^a?c 1 2 —a?c3 



The values for c n and therefore for a n can actually be 
determined as the Fourier coefficients of a certain uV(<£) 
curve (reference 3) . 

If the radius of the circle in the Z plane is denoted by 
R, tho coordinates of any point on this circle can be 
expressed as 

X=R cos £ 
Y=-R sin £ 

where, as £ increases from 0 to 2-tt, the" point describes the 
circle in a clockwise sense, thus leaving the external region 
to the left of the direction in which the boundary is trav- 
ersed. Any point of the circle can also be expressed as 

Z=X+iY=Re~ u 

If f =£+i77, the transformation 

Z=Re- f : 

yields the circle of radius R for 77=0 and, furthermore, 
7j= co corresponds to the region at infinity of the Z plane. 
Upon substitution of this expression for Z in equation (6), 
it follows that 

z=Re'c-*+c 1 +^e*+^-Wt+^e-*>e 3 *+ ... (9) 



10 



REPORT NO. 752 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



When this equation is separated into its real and imaginary 
parts, it is seen that 



x=Ci+Ee" cos £ +^e~* cos £+jpe~ 2 ' cos 2£ 



e" 3 * cos 3f + 



a=—Re" sin sin f +^c~ 2 ' sin 2f 



(10) 



+jge- 3 * sin 3{+ 



where the quantities Ci, a i} Oj, 03, . . . are deterrnined 
according to the method of reference 3. 

Rather than as a conformal transformation of a plane z 
into a plane Z, equations (10) are to be looked upon as the 
equations of transformation* from the rectangular Cartesian 
coordinates (z, o) into the orthogonal curvilinear coordinates 
(£, 17). Furthermore, the coordinate line ij=0 is the profile 
itself; that is, when -j=0, equations (10) yield the para- 
metric equations of the profile in the plane 2. It will he 
seen in the following discussion that the use of this type of 
coordinate system leads to a simple form for the boundary 
condition at the solid surface. 

BOUNDARY CONDITIONS 

If a body of revolution moves with a velocity U in the 
direction of its axis of symmetry, the normal velocities of the 
body and the fluid in contact with it are the same. Thus, 
at the surface of the moving body, the boundary condition 
expressed in the conformal orthogonal coordinates (£, -7) is 
simply 

since the coordinate lines along which 77 varies and £ remains ' 
constant are normal to the boundary 17=0. According to 
the first of equations (10) this boundary condition can be 
written as 

-QrXo =u ( R 008 *"3 003 £ 

cos 2£-^f cos 3|- . . .) (11) 

Furthermore, as the fluid at infinity, originally at rest, 
remains undisturbed by the motion of the body, the bound- 
ary condition there is, simply, 



(ML- 



(12) 



It is remarked that the simple form of the boundary con- 
dition (11) has been attained on account of the introduction 
of conformal orthogonal coordinates (£, ij) so defined that 
one of the coordinate lines 17=0 is the meridian profile itself. 
It will be seen in the following discussion that this choice of 



coordinate system leads to a process of iteration for the solu- 
tion of the fundamental differential equation (5) involving 
only simple quadratures. 

SOLUTION OF THE FUNDAMENTAL DD7FERENTIAL EQUATION 
BY ITERATION 

If the right-hand side of equation (10) for a is substituted 
into equation (5), it follows that 



(13) 

With regard to an iteration method it is desirable that the 
initial step in the process be obtained in a closed form. 
Thus, equation (13) has been so arranged that the solution 
of the left-hand side placed equal to zero can be obtained in 
a closed form. This initial solution can then be utilized as 
the starting point of an iteration process. Before a detailed 
description of the iteration method is given, however, it is 
first necessary to introduce several new parameters. Thus, 
the coefficient Oi is replaced by s 2 ^ and the radius B by ae" 
where a, as in equation (7), serves merely to preserve dimen- 
sions. Then equation (13) can be written as follows: 



003 *b| 



sinn {v+cc-fi sin ^^+^)+sinh fo+a-,3) 
+cosh (v+cc-p) sin ^=X>, e -c»+J)(-+a-fl) 

X[sin (n+l){(^+^) + (7H-l)cos(n+l)^ 
-( W +l)sin(n+l)f^] (14) 

where E»= 2 (<^g)»+2 > a ' nondimensional quantity, and tho 
boundary condition (11) becomes 

~r 2aUet> [sinl1 (a_/3) °° s e 

— 2V" 2(a_?) cos 2f— 36je- 3<a ~« cos 3f— . . . ] (16) 

The method of solving equation (14) is based upon the 
following considerations: 

It is assumed that the velocity potential <t> can be developed 
in a series of terms each of which is homogeneous with 
respect to the indices of the coefficients b n ; that is, 



4>-- 



n-0 



(16) 



where, for example, 
^0 involves none of the coefficients b„ 
4>i contains the coefficient 61 as a factor 
& is the sum of two terms having, respectively, bi 2 and 63 
as factors 

#8 is the sum of three terms having, respectively, &1 3 , b\bi, 
and 63 as factors and so on for the higher-order terms 
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Equation (16) for <f> is now substituted into equation (14) 
and the terms involving the coefficients b H and their products 
to the same degree in the indices are equated. The same 
process is applied to the boundary condition (IS). In this 
manner, equation (14) is replaced by a set of partial differ- 
ential equations which, in conjunction with the boundary 
conditions, can be solved rigorously for ft,, </> u <fo, . . . When 
the operations just described are performed, the differential 
equations with the accompanying boundary conditions for 
the first three functions <fo, ft, and <fe are as follows: 

sinh („+«-# sin ?(^+^)+sinh (,+a-/3) cos*^ 



+cosh(„+a-/3)8inf^=0 

with the boundary condition at the surface 



(17) 



and at infinity 

sinh (H-a-fl sin $ @+^) 
+sinh (ij+a— j8) cos f^+cosh (ij+a— /3) sin £^ 1= 

+2 cos 2^-2 sin 2^] 
with the boundary condition at the surface 

and at infinity 

sinh (,+«-« sin*(^+^) 



(18) 



+sinh (ij+a— jS) cos{-g|+cosh (ij+a— 0) sin$-=^= 



b&_ 



bu 



6 ie -^^)[sin 2* (|£ + |fc)+2 cos 2^-2 sin 2^ 1 ] 
+6 a e-3<v+-/>»[sin 3 f @+^)+3 cos 3^-3 sin 3«|£] 



(19) 



with the boundary condition at the surface 

(t^) 0 = 6aZ7e ^ e_3(a ~" ) cos 3£ 



and at infinity 

DETERMINATION OF <fo 

In order to solve equation (17) for <fo, it is convenient to 
introduce a new set of independent variables fi and X, where 

n=coa £ and X=cosh (ij+a— jS) 

By use of the symbolic expressions 



b . ,b 

2)* .b . . . . b* 

_ = _ C08 f _ +smS ^_ 



and 



^=sinh (u+a-P)^ 

b s b b* 

^5= cosh (ij+a— j8)-^+smh 2 fo+a-tf)^ 



bX* 



(20) 



b 

5i 



it follows that equation (17) can be written as 

i[t-^]-ih!]- w 

If asolutionof theform <fo=F(ji) G(\) is assumed, equation (21) 
becomes 

Since the left-hand side of this equation contains /t but not X 
and the right-hand side contains X but not p, F(ji) and #(X) 
must be such that each side is a constant. If, furthermore, 
the constant is chosan to be — 71(91+1), where n is a positive 
integer, then 



l£ (i -^i] +n(n+i)i M 
^[( i - x ^+-^+ i > Gr =°j 



(22) 



Each of these equations is a Legendre differential equation 
and, being of the second^ order, possesses two independent 
particular solutions: namely, the Legendre functions of the 
first and second kinds, P„ and Q„, respectively. According 
to the boundary condition at the surface, the form to be 
taken for P(/i) is Pi GO- Correspondingly, according to 
the boundary condition at infinity, the form to be taken for 
Q (X) is Qi(X), since the Q n (X) and their derivatives vanish 
for X= oo . It follows that the solution of equation (21) is 

<h=A 1 P 1 ( fi )Q 1 (\) (23) 

The arbitrary constant A\ is determined by the boundary 
condition at the surface: namely, 
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Thus, 



where 



and 



2aZ7e" 



<2i'(\>) 



X 0 =cosh (a—P) 

It will be necessary for the determination of fa and fa to 
have the first few orders of P K (ji) and <?„(X) stated explicitly. 
The following arc expressions for the first four harmonics: 

PoW = l 
Pi(m)=m=cos£ 

P 2 W=|(3M 2 -D=|cos2f+i 
' P s (m)=-|(V-3/x)=|cos3{+|cos{ 
QoW^log-^y 
Qx(X)=-ixlog^±-J-l 
Q a (X)=i(3X'-l) log^-|x 

<2,(x)=i(5X 3 -3x) iog|±Y-|x s +| 

DETERMINATION OF in. 

When the independent variables p. and X are introduced and 
equations (21) and (23) are used, it follows that equation (18) 
for fa can be written as 

£[^f]-!{ (1 -^]= 
|6 1 ^J^^^U2$ 1 (X)-(X+2VV=I)(2/(X)]P 0 (, l ) 

-2[Q,(X) + (X+2V>?=l)(2i'(X)]P 2 W} (25) 

The right-hand side of this equation requires that 

*i=Po(X)Po(p)+P 2 (X)P2(m) - (26) 

If this expression for fa is substituted into equation (25) and 
the coefficients of P 0 (n) and P 3 (ji) are equated on both sides 
of the resulting equation, it follows that 

|6 1 A 1 ( X -^7 I I 1 [(X+2- v /X^T) Q 1 '(X)-2Q l (X)J (27) 

6 yx — l 

^[(l-X^]+6P a (X) = 

|^ 1 ^ fr- / ^E ! -l 2 [(X+2VX^l) &'(X) + <2i(X)] (28) 



and 



With the right-hand sides taken equal to zero, these two 
equations reduce, respectively, to Legendre's equation for 
the zero- and second-order functions. Then, according to 
the theory of linear differential equations of the second order, 
the solutions of equations (27) and (28) are: 

Po(X)=P 0 &(X)-J^JPo(X)dX (29) 

and 

P S (X) =P 2 &(X) - (3X S -1) J^z^jz^p J (3X>- l)i2,(X)dX 

where P 0 (X) and P^X) are, respectively, the right-hand sides 
of equations (27) and (28). It is noted that each of these 
solutions contains only one arbitrary constant, B 0 and B 3) 
because the other two independent solutions Po(X) and 
P S (X), respectively, of the homogeneous forms of equations(27) 
and (28) do not satisfy the boundary condition at in- 
finity and therefore the arbitrary constants associated with 
them are taken equal to zero. 

The boundary conditions for fa lead to the following 
boundary conditions for P 0 (X) and P 2 (X): 



and 



(31) 



(32) 



The integrations required in equations (29) and (30) are 
straightforward and need not be performed hero. Further- 
more, by use of the well-known recurrence formulas for the 
Legendre functions, it can be shown that 

FoCK) = -§Mi {f [2<2o(X) + & (x)]+ yx 5 ^ [_Qo (X) - Qi (X)]| 

(33) 

and 

P 2 (x) =P 2 &(x) {^[2<2o(x) +&(x)_ 
+ yl?= r l[<2o , (x)+2Q 1 (x) 



(34) 



where, by means of the boundary conditions, equations (31) 
and (32), 



Po=0 



B 2 =66 1 ^ 1 ^^{(2 1 / (Xo) +|yV z: ^[«o , (x,)] 2 }] 



(36) 



It may be easily verified that these expressions for F 0 (X) 
and Fi (X) satisfy the condition that they and thoir deriva- 
tives vanish for X= 00 . 
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When the variables /i and. X are introduced into equations (21) and (25), it follows that equation (19) for fa can bo 
written as 

+i/-^^ , [8,(l-^-(V-l)(X + 3V^)t]+«. ! ^^^'[Kl-,^-/(X+ 2 Vx^l)|*] (30) 

If the expressions for fa and fa given, respectively, by equations (23) and (26) are substituted in equation (36), it 
follows that: 



(37) 



where 



and 



B 1 (X) = -26 I ^^^^^ 2 (X)-(X+2Vx^l)[ J F'o / (X)+^'(X)]) " 

+ |6Mi (X ~^ l) ^ 1 (X) + (X+2Vx^l)^(X)] 

The right-hand side of equation (37) requires that 

fa=F 1 (\)P 1 Qi)+F 3 (\)P i ( fl ) (38) 

When this expression for fa is substituted into equation (37) 
and the coefficients of Pi Qi) and P 3 (X) are equated on both 
sides of the resulting equation, the following equations for 
F x (X) and F 3 (X) are obtained: 

k[G-*)^]+WiW^*.(X) (39) 



and 



A VJ^J+^W^OO (40) 

It is evident that the homogeneous forms of these equations 
are, respectively, Legendre's equations for the first- and 
third-order functions. According, then, to the theory of 
linear differential equations of the second order, the solutions 
of equations (39) and (40) are 

F l (.\)=B 1 Q 1 (\)-J^j^- T) jR 1 (\)d\ (41) 



and 



F 3 (\)=B 3 Q 3 (\) 

-(5X 3 -3X)J (x2 _ 1) ^ X xS _ 3X)2 J(5X 3 -3X)ii! 3 (X)dX (42) 

The arbitrary constants B t and B 3 are determined by the 
boundary condition at the surface of the moving body. 
Thus, the boundary condition for fa leads to the following 
conditions for Fi (X) and F 3 (X) : 



and 



(43> 



(44) 



Then, if the integrations required in equations (41) and (42) are made and the recurrence formulas for the Legendre 
functions are used, it can be shown that 

Fx (X) =B 1 Q 1 (X) -| 6 A j| Qi (X) + 1 & (X) + Vx^IQ/ (X) - Q t (X)]J 

- 1 bMi j-4 -|§ Q l (X) -| &(X) +| Q 0 ' (X) + - V / X 5 =l [Q t (X) +4Qo(X) - 5&' (X)]J 

-| Mi jl + J «i(X) +f &(X)- VJ?=T [&(X) + Q 0 (X) -| (X)]} (45) 
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F s (\) =5 s ft (X) -§ b& (| ft (X) +| ft (X) + Vx^fft'W +4ft (X)]] 

- 1 -§ ft (X) +i£ ft (X) +\ ft' (X) + V?? 11 ! [- 6ft (X) + ft (X) -§ ft ' (x)] ) 

|+| ft (x) ft(x)+ V>?=i [§ft(x) -| ft(x)+ft'(x)] j 



8 A d 



(46) 



where 



and 



Sift' (Xo) =| 6x5, (3ft (Xo) + ft'{Xo) -| -v^T (3ft (Xo) + ft' (Xo) -§ [ft' (Xo) ]' }) 

-| & i 2 ^(3ft (Xo) +±± ft' (Xo) -f ft" (Xo) -§ VV^T {3ft (Xo) + ft' (Xo) [ft' (Xo)] 2 }) 
-|6 2 ^ 1 (-3ft(X)+^-ft'(X a )-VV Z: l {3ft'(Xo)-i[ft'(X 0 )] s J) 

-B s ft' (X,) =| b& (3ft (Xo) + ft' (Xo) --^J==[^ ft (Xo) ft (Xo) + ft' (Xo)] } 

§ v.*, f£ ft ex,) +| ft'(Xo) +| ft"(Xo) -y=rr [-ft & o*> +f & w -f &'(*a>]) 



(47) 



+ 



+| Mi ft (Xo) +^ ft' (X,) -8ft(Xo) --WT^J t&' -2]) 



(48) 



It may be easily verified that the expressions for Ji(X) 
and J's (X) and their derivatives vanish for X= « . 

It is clear that the iteration can be continued to obtain 
further members of the sequence of The order n to 
which the process must be carried depends on the magnitude 
of the coefficients b n of the conformal transformation. For 
example, it may be that one of the higher coefficients, say 
b a , is significant. It is obvious, then, that the iteration must 
be continued far enough to include at least the corresponding 
term <j> m . Although the <£ B 's do not form a set of normal func- 
tions (obtained independently of one another), they are 
nevertheless uniquely determined by means of simple quad- 
ratures; that is, they constitute a set of universal functions in 
the sense that the labor necessary for their determination is 
performed but once and need not be repeated for each, par- 
ticular body of revolution. 

VELOCITY AND PRESSURE 

According to Bernoulli's theorem, 



where 



(49) 



p pressure anywhere in fluid 

p* pressure in region at infinity undisturbed by motion of 
solid 

p density of fluid 

g r velocity of fluid relative to boundary of solid 



It is recalled that in this paper the expression for the 
velocity potential <f> has been derived by considering the 
coordinate axes as rigidly attached to the body and there- 
fore moving with a velocity U in the positive direction of 
the x axis. The velocity vector 2= — grad <j> represents the 
velocity of the fluid relative to the undisturbed fluid at 
infinity. The vector velocity q_ T , relative to the moving 
axes, is therefore given by 



or 



q T =i- U 
2r=t(2z-C7)+J2« 



where i and j are, respectively, unit vectors along the posi- 
tive directions of the x and a axes. The magnitude of the 
relative velocity is therefore obtained from 



where 



2*= — || and g, 



(50) 



da 



It is recalled, however, that the velocity potential <f> is 
derived as a function of the independent variables £ and rj. 
In order then to determine g r , it is necessary to express the 
velocity components q x and in terms of the quantities 



d<f> , d£ 
and 5— 

lowing way: 



These expressions can be derived in the fol- 
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The relations 



z=x+i© and z=x — io 



can be considered as equations of transformation from the 
coordinates (x, o) to the independent coordinates (z, z). It 
follows then, symbolically, that 



and 



Therefore 



2 a =A_^ 

bz dx ba 



2 dz dx^do 



2x+t 2o =-(| r H^)*=-2!! 

where <f> is now a function of the independent variables z 
and 2. 

The conformal transformation z=/G") then gives 

. . n b<i>df 



or 



where 



Symbolically, 



and therefore 



, . 2 dz b4> 
2*+^ J*Ttdi 

_ dzdz 



n cZz /dx . d©\ . . /do dx\ 
Then, by means of the Cauchy-Riemann relations 



it follows that 



or 



dx &o , ite bo 

-sr^ ^ bi, — &i 

dz_bx .bx 
<Zz_ do . .So 



(51) 



Equation (51) therefore becomes 

, . 1 /bx .bx\/'b<j> , .d<A 

— HT^-dy Uf+^J (52) 

where 

MINIKINS 



When the real and imaginary parts on both sides of equation (52) 
are equated, it follows that 



1 /d0 dx d# bx\ 

1 /b<t> ba Oo\ 

2o_ J s Vdf d?" 1 "^ a,/ 



(53) 



By means of these equations q/ and therefore the pressure 
coefficient f — — given, respectively, by equations (50) and (49) 

can be obtained as functions of the independent variables 
£ and ij. 

APPLICATION OF THEORETICAL RESULTS 

SYMMETRICAL JOUKOWSKI SHAPE 

General expressions. — The theoretical results of this paper 
will now be applied to the case of a body of revolution whose 
meridian section is a symmetrical Joukowski profile with- 
both rounded nose and tail. This example is sufficiently 
complicated to illustrate the principles and usefulness of 
the method. 

It is well known that by means of the mapping function 
z=Z'+^r (54) 

the circle of radius a, with its center at the origin of the Z' 
plane, is transformed into the line segment extending from 
z=— 2a to z=2a in the z plane; and the circle of radius 
a(l+«i+«s) with its center at Z'=t& is transformed into a 
symmetrical Joukowski profile with rounded nose and tail 
in the z plane (fig. 2). When Z' is replaced by Z+eja, 
equation (54) becomes 
TY 



Ada 






iT o 







-XX 



8a,t 




2a,C 



Z'.Z-phnes- z-plane 
Fiouee SL— Mapping of a circle Into a symmetrical Joukowski profllo. 



z=Z+e l a- 



or 



Z=Z+€iffi-f-^— €i^+£! S ^— «l 3 ^i+ • • • (55) 

A comparison of this equation with equation (6) shows that 
c l =e 1 a, ai=a 2 , a 2 =—6 1 a? l a s =e 1 1 a i , at—— efa 6 , . . . 

Also, the radius of the circle in the Z plane being denoted 
by ae", it follows that 

««=l + erf-* (56) 



Z+^a 
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According to definition 



where 



As a t =a 2 and 0*+!= (— l) ra e 1 e, a. m+2 , it follows therefore that 



/J=0 and &„=2(-l) n ei" 

The expression for the velocity potential then becomes 

<£=<fo+<£i+<fe+ • • • 
where, from equation (23) 

with 



(57) 

(58) 
(59) 



A 1 =- 



2aU 



ft'(Xo) 



and 

Xo=cosh tt=s |(i +6l +^+ 1+ ^ + J 
From equation (26) 

fc=JP,(A)P,(lO+*«<A)PiG«) (60) 

where 

P„(X) (§[2ft (X) + ft (X)]+ Vx^ltft' (X) - ft (X)]J 

P 2 (X)=.B 2 ft(X) 



+| ^ J [2ft (X) + ft (X)] + V??=l [ft' (X) +2ft (X)] 



and 



Xo 



ft'(Xo) 



3ft , (X 0 )+2VV i -l[ft / (Xo)] 2 } 



From equation (38), 
where 



and 



0 s =P 1 (X)P I a t )+P 3 (X)P 3 (M) 
Pi(X) =Pxft(X) +| eA {§ Pft (X) +3ft(X)]+ VX^T [ft' (X) - ft (X)] J 

+| e,% { 2+| ft (X) -| ft(X) -| ft'(X) + VX*=I [~2ft(X) + ft (X) +^ ft' (X)] ) 

P 3 (X)=B,ft(X)+| ei B 2 [I [2ft (x)+3ft(X)]+ Vx*=T [ft' (X) +4ft(X)] j 

{-£-§ ft (X) -i| ft(X) -| ft'(X) + V>?=t[| ft(X) + fft(X) +\ ft'(X)] j 

Bift'(Xo) = -f «A (3ft(X 0 ) + ft'(Xo) -| VxF : l{3ft(X 0 )+ft'(X,) -f [ft'(Xo)] 2 }) 

(-6ft (Xo) +3ft (Xo) +^ ft' (Xo) - J ft"(Xo) — VX7="1 { 9ft (Xo) + 15ft'(Xo)-ll[ft'(Xo)] 2 }) 

5 3 ft' (Xo) = -f 6,5, { 3ft(Xo)+ ft'(Xo) -j=T=i & W ft (X,) + ft' (X,)] ) 

-fftA, (l6ft(Xo)-^ft(Xo) -^i'(Xo) -ift"(Xo) - [l6+f|ft (Xo) -^ft(Xo) -^ft'(Xo)] ) 



(61) 



Velocity and pressure. — From equation (54) together with 
the transformations 

Z'=£ia+Z and Z=ae a e~^ 



it follows that 

cos {) (i+ ^ +2 ^.^ s g+?gra ) 



q=— ae*" 1 "" sin £ 
Then 



0-5 



5 +2l^e^ T;; ^o¥l+e 5 <^ 



) 



(62) 



&{~ <wTfttsm fj 1 [^+2^" cos S+e 2 '^] 2 ) 



and 



, [e^+e^^e^" cos {+2f 1 e 3f " + '' ) 

&, -a « acosf [ ei 2 +2 ei e<+« cos S+e 2 '"™] 2 



At the surface of the body i?=0, so that 



. fa'+ g*'")*' cos g+2 6l e a 



■] 



(63) 



Therefore 



(J*)_o=aV«[l- 



l+2e t °(l-cos 2£) 



"£ 1 2 +26 1 e a cos £+e 2an (e 1 2 +2eie a cos f+e 2 ' 



(64) 
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&(Xo) =0.82657; &W=- 15.99664 
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The boundary condition at the surface yields immediately 
the expression for o ' Thus 



5 =o 



(65) 



where the expression for (^f ) 0 ^ £i veu Dy tae second of 
equations (63). 

The expression for (j^J is obtained from the equation 
for <£. Thus, the velocity potential can be written as, 



*=4,P, GO Q t (X) + [F 0 {\) P„G0 +F,(\) P.O.)] 
+ [P 1 (X)P 1 0x) + P 3 (X)P3(m)]+ • • • 



Then 



-(||) i _ o =[4 l «i(X a )+P,(Xo)+|P3(Xo)]sin f +|P 2 (Xo) sin 2? 
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+^P 3 (Xo) sin 3|- 



(66) 



where A x , Pi(Xo), Pi(Xo), and P 3 (Xo) are obtained from 
equations (59), (60), and (61). 
By means of equations (63), (64), (65), and (66) the 

velocity g> of the fluid and the pressure coefficient -j 12 

\ P TP 

given, respectively, by equations (50) and (49) can be 
evaluated. 

Numerical example. — As a numerical example, consider the 
case where «,=0.15 and «j=0.10. Then e a =l+( 1 +e 2 = : -^ 
and 



/x_\ 209 5 
V2oA-o~120 + 8 C01 



1540 



1 + 



3 317+75 cos £ 
200 



Also 



317+75 cos £ 



61600 



) 



(67) 



(317+75 cos £) 



5(IL""8 ai,lt [ 1 - 

Wd£\ 5[~ 200( 75+317 cos £) ~| 

2a\c)7,/,.n == 8L COS 5 (317+75 cos £)' J 



(68) 



and 



\2a /»-o 64L 
Now, 



400 



317+75 cos £ 



5000(33-25 cos 2 £)"| . 
(317+75 cos £)* J ^ 



Xo= — 2 — = 40 VV— 1 = 



e g — e" 
2 



9_ 
= 40 



Therefore, from equations (24) 
<2o(Xa) =2.19723; &'(\») = 



1600 • 0'V>^- 5248 
-§r» «o W-^561 



&(Xo) =0.57729; &'(Xo) = -13.91684 
and from equations (59), (60), and (61) 

Ai = 0.055403 (2a U) 

P, = 0.502158e 1 i (2at7) 

B 2 = 0.43 1087ei (2a U) 

53=0.537892^(2(1 E7) 
P 1 (X 0 ) = 0.079143e 1 2 (2aC7) 
P 2 (Xo) =0.2451 19ei(2aZ7) 
i? 3 ( =- o. 180438« 1 2 (2o V) 

Then, from equation (66), 

0 ^ 2aC7 I(°- 069373 + 0 - 011479ei2 ) sinf 
+0.367679 £l sin 2£-0.338321 £l s sin 3£] 

or with e 1 = 0.15, 

o =2aZ7(0.069631 sin £+0.055152 sin 2£ 
"°- 0.007612 sin 3£) 



From equations (65) and (68), 



200(75+317 cos £) 
(317+75 cos £)' 



] 



(70) 



(71) 



Table I presents the calculated values of the quantities 
given by equations (68), (69), (70), and (71) for intervals of 
10° in the angle £. Table LT gives the coordinates of the 
profile and the corresponding values of the velocity q T /U 

and the pressure coefficient K — — • Figure 3 shows the 




Fiouee 3. — Pressure distribution at the surface ol a JoukowskI shape. 

graphs of the meridian profile and the pressure distribution. 

The small circles on this graph denote the values of 

\oTT- 

obtained by the method of reference 2. It is seen that the 
numerical results obtained by the two methods agree almost 
perfectly. 
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LOW-DRAG SYMMETRICAL AIRFOIL SHAPE 

The abscissas and ordinates of the meridian profile are 
given in table HI. The first step is to obtain the^(<£) curve 
of the meridian profile according to the method of Theodor- 
sen and Garrick (reference 3). The <f> in reference 3 corres- 
ponds to the — £ of the present report. A Fourier analysis 
of the^(<£) curve is then performed and immediately yields 
the conformal transformation of the meridian profile into a 
circle. This conformal transformation is then looked upon 
as the equations of transformation from the rectangular 
Cartesian coordinates (x,a) of the meridian plane into the 
orthogonal curvilinear coordinates (£, 17), denned in such a 
manner that the coordinate line ij=0 is the meridian profile 
itself. For t;=0, then, these equations yield the parametric 
equations of the profile and for the present numerical 
example they are as follows: 



^=0.08895+0.98698 cos £— 0.05586 cos 2£ 
+0.03553 cos_3£— 0.01122 cos 4£+ . . . 

|j=- 0.25762 sin £- 0.05586 sin 2£ 

+0.03553 sin 3£-0.01122 sin 4£+ . . . 

The chord c of the profile is given by 

c=2.0450a 



(72) 




O 10 20 30 40 50 60 70 80 90 100 
Percerff chord 

Fiqdbk 4.— Pressure-distribution curves. 



TABLE I 

CALCULATED VALUES FROM EQUATIONS (68,) (69), (70), 

AND (71) 



(dig) 


20(51),-° 


2a(e*),-o 








0 


0 


a 30612 


0 


-a 30612 


0.09371 


6 


-.07636 


.30426 


-.01368 


-.30426 


.09840 


"10 


-. 16220 


.29871 


-.02715 


-.29871 


.11242 


16 


-.22746 


.28963 


-.04022 


-.2S953 


.13667 


20 


-.30146 


.27680 


-.0S267 


—.27688 


.16763 


30 


-.44446 


.24178 


-.07497 


-.24178 


.25600 


40 


-.67824 


.19543 


—.09248 


-.19543 


.37266 


GO 


-.69990 


.14049 


-.10385 


-.14049 


.60959 


60 


-.80658 


.0S02J 


-.10807 


-.08025 


.66702 


70 


-.89544 


.01849 


-.10469 


-.01849 


.80216 


SO 


-.96362 


-.04072 


-.09403 


.04072 


.93022 


90 


— L 00813 


-.09329 


-.07724 


.09329 


1.02502 


100 


-L 02684 


-.13662 


—.05630 


.13552 


1.07071 


110 


-1.01361 


-. 16466 


—.03379 


.16455 


1.05428 


120 


-.96807 


-.178S9 


-.01254 


.178S9 


.96917 


130 


-.88699 


-.17896 


.00478 


.17896 


.81876 


HO 


-.76911 


-. 16736 


.01615 


.16735 


.61953 


ISO 


—.61662 


-.14856 


.02056 


.14866 


.40096 


1C0 


—.43043 


-.12888 


.01823 


.12888 


.20188 


170 


—.22162 


-.11399 


.01058 


.11399 


.06211 


176 


-.11164 


-.10988 


.00548 


.109SS 


.02454 


ISO 


0 


-.10847 


0 


.10847 


.01177 



TABLE II 

RELATIVE VELOCITY AND PRESSURE COEFFICIENTS 
CORRESPONDING TO PROFILE COORDINATES 



(dog) 


X 

2o 


a 

2a 


1' 
t? 


p— Pa 
IP V 


0 


1.0572 


0 


0 


1 


5 


L 0638 


. 0267 


.28700 


. 91766 


10 


1.0438 


. 0530 


. 63619 


.71357 


15 


L0273 


0787 


.72700 


' 47148 


20 


L0042 


!l034 


! 88521 


.'25141 


30 


.9390 


.1489 


1.02660 


-.05390 


40 


.8496 


.1872 


1.09887 


-.20761 


50 


.7378 


.2165 


1.12592 


-.28769 


60 


.6061 


.2369 


1. 12840 


—.27330 


70 


.4673 


.2446 


1. 11667 


-.24606 


SO 


.2948 


.2425 


1.09660 


-.20253 


90 


.1223 


.2307 


1.07204 


-.14927 


100 


-.0556 


.2105 


1. 04680 


-.09309 


110 


-.2340 


.1842 


1.02000 


-.04030 


120 


-.4074 


.1640 


.99609 


.00781 


130 


-.6699 


.1225 


.97600 


.04933 


140 


-.7149 


.0922 


.96652 


.08488 


150 


—.8383 


.0546 


.93906 


. 11718 


160 


-.9279 


.0403 


.91741 


.16830 


170 


-.9861 


.0193 


.84684 


.28287 


176 


-.0997 


.0096 


.67771 


.64071 


ISO 


-L0046 


0 


0 


1 



TABLE HI 

ABSCISSAS AND ORDINATES OF MERIDIAN PROFILE 



I 


a 


(percent chord) 


(percent chord) 


0 


0 


.60 


XOS 


.75 


2.53 


1.25 


3.25 


2.50 


4.62 


6.00 


6.69 


7.50 


8.06 


10.00 


9.23 


16.00 


11.04 


20.00' 


12.32 


26.00 


13.23 


30.00 


13.84 


36.00 


14.20 


40.00 


14.29 


45.00 


14.13 


50.00 


13.62 


65.00 


12.69 


60.00 


11.34 


65.00 


9.78 


70.00 


8.03 


75.00 


8.24 


80.00 


4.47 


85.00 


2.88 


90.00 


1.66 


95.00 


.54 


100.00 


0 



Equations (72) may then be written as 



|=0.04350+0.48263 cos £-0.02732 cos 2£ 
+0.01737 cos 3£- 0.00549 cos 4£+ . . . 

®=— 0.12598 sin £—0.02732 sin 2£ 
+0.01737 sin 3£-0.00549 sin 4£+ . . . 



(73) 



Table IV gives values of — and — calculated according to 

c c 

these equations and figure 4 shows the very good agreement 
between these points (small circles) and the profile plotted 
from table LIT (solid line). 

"With the coefficients of the conformal transformation 
known, the next step is to obtain an expression for the veloc- 
ity potential <f> at the surface of the body. For the present 
numerical case, 



<2>=2aZ7(0.00687+ 0.08181 cos £+0.03384 cos 2£ 
— 0.02251 cos 3£+ . . .) 



(74) 
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TABLE IV 

CALCULATED VALUES PROM EQUATIONS (72) AND (73) 



( 


tint 


«In2f 


sln3f 


sin 4? 


cos { 


COS2f 


cos 3£ 


COS4f 


x 

— 
c 


oj 

C 


Q 


U 


Q 


o 


— 

o 


1 


1 


1 


1 


0.5107 


0 




08716 


. 17385 


.25382 


.34202 


.99619 


.98481 


.98693 


.93969 


.6090 


— . 0131 


10 


\ 17386 


.34202 


Isoooo 


.61279 


1 98481 


193969 


.86603 


.76604 


.6010 


—.0261 




. 25882 


. 60000 


. 70711. 


1 86603 


196593 


.86603 


.70711 


.50000 


.4956 


—.0388 




34202 


. 64279 


.86603 


198481 


1 93909 


1 76604 


.60000 


.17366 


.4838 


—.0610 




. 60000 


. 86503 


1.0000 


\ 88603 


. S6603 


.'5OOOO 


0 


—.50000 


.4606 


—.0740 


40 


64270 


198481 


' 86603 


*34202 


' 76604 


1 17365 


—.60000 


—.93969 


.4049 


—.0947 


fiO 


;76601 
.86603 


198481 


Ifioooo 


-.34202 


.64279 


-.17385 


-.88603 


-.93969 


.3486 


-. 1129 


00 


.88603 


0 


—.86603 


.50000 


-.60000 


-1.0000 


-.60000 


.2838 


-.1280 


70 


.03960 


.64279 


-.50000 


—.98481 


.34202 


—.76604 


—.86603 


.17366 


.2135 


—.1392 


80 


.98481 ■ 


.34202 


-.88603 


-.64279 


.17385 


-.93969 


-.50000 


.76604 


.1401 


-.1449 


90 


1.0000 


0 


-1.0000 


0 


0 


-1.0000 


0 


1.0000 


.0653 


-.1434 


100 


.98481 
.93969 


- -.34202 


-.88603 


.64279 


-.17365 


-.93969 


.60000 


.76604 


-.0102 


-.1333 


110 


. -.64279 


-.60000 


.9S4S1 


-.34202 


—.78604 


.88603 


.17356 


-.0866 


—.1149 


120 


.86603 


-.88603 


0 


.86803 


—.50000 


-.60000 


1 


-.60000 


-.1640 


-.0902 


130 


.76604 


-.68481 


.60000 


.34202 


—.64279 


-. 17365 


.88603 


-.93969 


-.2418 


-.0628 


140 


.64279 


-.98181 


.88603 


-.34202 


-.76804 


.17386 


.50000 


-.93969 


-.3171 


-.0372 


1£0 


.60000 


-.88603 


1.0000 


-.88603 


-.86603 


.50000 


0 


-.60000 


-.3864 


-.0172 


160 


.34202 


-.64279 


.86003 


-.98481 


-.93969 


.76604 


-.50000 


.17365 


"-.4406 


-.0051 


170 


.17386 


-.34202 


.50000 


-.64279 


-.98481 


.93969 


-.86603 


.76604 


-.4787 


-.0003 


ISO 


0 


0 


0 


0 


-1.0000 


1.0000 


-1.0000 


L0000 


-.4893 


0 



TABLE V 

CALCULATED VALUES FROM EQUATIONS (49) TO (53) 





1/- 


S/* 


(£)' 








id 

U 


ft)' 






0 


0 


0.3077 


a 0947 


0 


0.3077 


1 


0 


0 


0 


1 


16 


-.2381 


.2927 


.1414 


.0073 


.2927 


.6937 


-. 6037 


.4188 


.6472 


.5812 


20 


-.3138 


.2822 


.1781 


.0130 


.2822 


.4241 


-.5178 


.6998 


.7745 


.4002 


30 


-.4644 


.2686 


.2825 


.0320 


.2665 


.1809 


-.4624 


.8755 


.9357 


.1246 


40 


-.6013 


.2279 


.4136 


.0608 


.2279 


.0372 


-. 3648 


1.060 


1.030 


-.0601 


50 
60 


-.7147 


.1963 


.6493 


.0956 


.1963 


-.0542 


-.2896 


1.196 


1.093 


-.1952 


-.7969 


.1671 


.6597 


.1295 


.1571 


-.1190 


-.2206 


1.301 


1.141 


-.3008 


70 


-.8466 


.1026 


.7272 


.1641 


.1026 


-.1660 


-.1412 


1.377 


1.174 


—.1771 


80 


-.8703 


.0274 


.7582 


.1622 


.0274 


-.1862 


-.0374 


1.406 


1.186 


-.4060 


90 


-.8804 


-.0668 


.7796 


.1493 


—.0668 


-.1629 


.0883 


L380 


1.166 


-.8602 


100 


-.8890 


-.1686 


.8188 


.1159 


-.1686 


-.0911 


.2069 


1.233 


Llll " 


-.2333 


110 


-.9018 


-.2682 


.8798 


. .0671 


-.2582 


.0070 


.2843 


1.067 


LOSS 


—.0670 


120 


-.9128 


-.3137 


.9313 


.0122 


-.3137 


.0937 


.3115 


.9185 


.9584 


.0816 


130 


-.9040 


-. 3195 


.9193 


-.0377 


-.3196 


.1481 


.3011 


.8163 


.9034 


.1837 


140 


-.8521 


-.2734 


.8008 


-.0726 


-.2734 


.1706 


.2662 


.7589 


.8709 


.2412 


ISO 


-.7357 


-.1897 


.6772 


-.0862 


-.1897 


.1710 


.2188 


.7330 


.8561 


.2670 


160 


-.5469 


-.0954 


.3071 


-.0740 


-.0954 


.1612 


.1466 


.7251 


.8515 


.2749 


170 


-.2917 


-.0220 


.0866 


-.0427 


-.0220 


.1612 


.0641 


.7248 


.8512 


.2761 


180 


0 


.0066 


.0000 


0 


.0056 


I 


0 


0 


0 


1 



The boundary condition at the surface of the moving body 
is simply 

Since ^' =— i fc follows from equations (72) that 

-Qfy o =2af7(0.25762 cos £+0.11172 cos 2£ 
-0.10658 cos 3£+0.04489 cos 4£- . . .) 
and from equation (74) that 

-Qfy fl =2aZ7(0.08181 sin £+0.06768 sin 2£ 



(75) 



—0.06752 sin 3£+ . . 
Also, from equations (72), 



•) 



(76) 



^=^=-0.98698 sin £+0.11172 sin 2£ 
—0.10658 sin 3£+0.04489 sin 4£— . . . 

^=-•^=0.25762 cos £+0.11172 cos 2£ 
—0.10658 cos 3£+0.04489 cos 4£— ... 



(77) 



Equations (75), (76), and (77) suffice to determine the pres- 
sure distribution over the surface of the body with the use 
of equations (49) to (53) contained in the section entitled 
"Velocity and Pressure." Table V presents the calculated 
values of the various quantities given by equations (49) to 
(53), the lost column containing values of the pressure 

coefficient V~'P'° - Figure 4 shows the graphs of the pressure 

distribution for both the body of revolution (solid curve) and 
the airfoil of infinite span having the meridian profile as 
cross section (broken curve). ^ 
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